o 

(N 

o 



Metal adatoms on graphene and hexagonal boron nitride: 
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Periodically corrugated epitsixial graphene and hexagonal boron nitride (^-BN) on metallic sub- 
strates are considered as perspective templates for the self-assembly of nanoparticles arrays. By 
using first-principles calculations, we determine binding energies and diffusion activation barriers of 
metal adatoms on graphene and /i-BN. The observed chemical trends can be understood in terms of 
the interplay between charge transfer and covalent bonding involving the adatom d electrons. We 
further investigate the electronic effects of the metallic substrate and find that periodically corru- 
gated templates based on graphene in combination with strong interactions at the metal/graphene 
interface are the most suitable for the self-assembly of highly regular nanoparticle arrays. 

PACS numbers: 36.40.Sx. 61.48.De, 68. 43. Be, 68.65. Cd 
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I. INTRODUCTION 

Two-dimensional graphene and hexagonal boron ni- 
tride (/i-BN) are attracting considerable attention due 
to their extraordinary physical properties and perspec- 
tive technological applications.-'- Epitaxial single lay- 
ers of graphene and /i-BN can be grown via the chem- 
ical vapor deposition route on a large variety of metallic 
substrates;^ In many cases, the epitaxial layers of these 
two materials are extremely well-ordered but reveal long- 
wavelength periodic corrugations, or Moire patterns, re- 
sulting from the lattice constant mismatchi^^— Such su- 
perlattices are considered as promising templates for the 
chemical self-assembly of periodic arrays of nanoparticles, 
with perspective applications in ultra-high density infor- 
mation storage, catalysis, sensing, etc.* The potential of 
this approach has already been affirmed by the successful 
production of regular arrays of nanoclusters with narrow 
size distributions!^"— 

One attractive property of this approach to nanoscale 
self-assembly is the high degree of customization. Indeed, 
the space of allowed chemical compositions involves three 
degrees of freedom: (i) deposited metallic nanoparticles 
formed by virtually any metal from the periodic table; 
(ii) a monolayer of either graphene or /i-BN; (in) metallic 
surfaces, one of the many (111) fee or (0001) hep surfaces 
of forth or fifth row late transition metals. Understanding 
the roles these three factors play in the self-assembly pro- 
cess is of paramount importance for the rational design of 
nanoparticle arrays with novel properties and functions. 

In this work, we aim at understanding the chemical 
trends in the binding and diffusion of individual metal 
adatoms, the initial step which largely predetermines the 
overall self-assembly process. In particular, by using first- 
principles calculations we perform a systematic study of 
the adatom binding and diffusion as well as of the elec- 
tronic properties of the metal adatoms upon deposition. 
First, we focus on the chemical trends observed for metal 
adatoms varying horizontally and vertically in the peri- 
odic table, when absorbed on free-standing graphene and 



/i-BN. Then, we model the electronic effects of periodi- 
cally corrugated epitaxial layers deposited on two metal- 
lic substrates representing the limiting cases of strong 
and weak interactions at the interface. 

The present paper is organized as follows. In Sec. 
im we describe our first-principles methodology and the 
adopted models. Section IIIII discusses the binding and 
the diffusion of adatoms on free-standing graphene and 
/i-BN. Section [IV] is devoted to the investigation of sub- 
strate effects. The conclusions are drawn in Sec.lVl 



II. METHODS 

In our approach, the electronic structure is described 
through the use of the Perdew-Burke-Ernzerhof (PBE) 
exchange-correlation density functional within density 
functional theoryf^^ The ultrasoft pseudopotentialsi^ 
used in the present study treat the d-electrons of all 
transition metals as valence electrons. The semi-core sp- 
states were also treated explicitly in the case of early 
third row elements (K-Fe). The one-electron valence 
wave functions and the electron density were described 
by plane-wave basis sets with kinetic energy cutoffs of 
30 Ry and 300 Ry, respectively^^ All calculations were 
performed using the spin-unrestricted formalism. We 
used the pwscf plane-wave pseudopotential code of the 
QUANTUM-ESPRESSO distribution^^ The convergence of 
the results with respect to the simulation parameters was 
systematically verified. 

The constructed models are based on periodic two- 
dimensional 3x3 supercells in combination with 4x4 k- 
point meshes. Potential energy surfaces for the free- 
standing graphene and /i-BN systems were studied by 
relaxing all atomic positions with the initial coordinates 
of the adatom set to the high-symmetry positions shown 
in Fig. [1] This procedure allowed us to determine both 
the potential energy surface minima and the transition 
state configurations connecting the neighboring equiva- 
lent minima, since symmetry was retained during the 



relaxation. The validity of this approach was verified 
through nudged elastic band calculations J^ For specific 
cases, we could compare our results with those of other 
investigations,^^'-- finding good agreement. 

We investigated the effects of local electronic-structure 
variations due to the presence of metallic substrates by 
considering the case of Co adatoms on epitaxial graphene 
and ft,-BN supported by lattice-matched Ni(lll) and 
Cu(lll) substrates. These metallic substrates corre- 
spond to regimes of strong and weak monolayer-substrate 
interactions, respectively.^*^ Under realistic conditions 
these metals do not produce Moire pattern due to the 
small lattice mismatchj^i^i"— However, through models 
involving such structures, we could study electronic ef- 
fects induced by the Moire superlattice using sufficiently 
small simulations cells and without introducing large lat- 
eral strains. A very similar methodology was also used 
in Ref. i4!. The considered model calculations consist of 
3x3 two-dimensional slabs composed of 4 atomic planes 
of metal with both surfaces covered by graphene or h- 
BN in order to avoid spurious dipole-dipole interactions. 
The slab configurations were fully optimized, including 
the distance between the metal layers<^ 



III. ADATOMS ON SUSPENDED GRAPHENE 
AND h-BN 




FIG. 1: (Color online) Potential energy surfaces of the Co 
adatom on freestanding (a) graphene and (b) /i-BN obtained 
by constraining the in-plane position of the Co adatom on a 
6x6 points mesh in the unit cell. The same energy scale is 
used in both plots and the energy is referred with respect to 
that of an isolated Co atom. The high-symmetry positions 
are indicated with labels. The arrows show the lowest en- 
ergy diffusion pathways connecting the local minima, h sites, 
through the transition state configurations, the b sites and ie 
sites in the case of graphene and /i-BN, respectively. 



First, let us compare the potential energy surfaces 
(PES) for a representative case of a Co adatom on free- 
standing graphene and /i-BN (Fig.[T|). For both monolay- 
ers the PES minima correspond to the hollow (h) sites. 
Interestingly, in the case of Co as well as for most of 
the other metal adatoms the distortions of the graphene 
and /i-BN lattices are very weak. Although the two 
PES involving the Co adatom are qualitatively similar, 
the binding energy on graphene are substantially higher 
(1.60 eV vs. 1.03 eV). The stronger binding to graphene 
is a systematic feature which can be understood consider- 
ing that the chemisorption of metal adatoms is governed 
by the following attractive contributions: (i) strong cova- 
lent bonding, and (ii) moderate electrostatic interaction 
as a result of charge transfer between the adatom and the 
monolayer. While graphene is a semimetal, /i-BN is an 
insulator with a band gap of w6 eVi^ The availability 
of low-energy electronic states in graphene results in a 
more efficient interaction with the adatom states and al- 
lows for charge transfer between the graphene layer and 
the adatom, thus introducing an electrostatic component. 
The preferential binding to the h sites is very common, al- 
though we find a number of exceptions, especially among 
the heavy transition metals. In particular, for Pd, Ir and 
Pt, the bridge (6) site is the lowest energy position on 
graphene. The nitrogen on-top sites (^n) correspond to 
the PES minima on h-BN for V, Ni, Pd, Ir and Pt. 

The lowest energy pathways connecting the neighbor- 
ing local minima involve transition states at b and Ib 
sites for the Co adatom on graphene and /i-BN, respec- 



tively. The diffusion activation barriers are again higher 
for graphene, 0.40 eV vs. 0.13 eV for /i-BN. However, this 
trend is not systematic across the periodic table and very 
often the diffusion on /i-BN is characterized by higher ac- 
tivation barriers [cf. Fig.[51Ja,b)]. Notably, the calculated 
activation barrier in the case of ft,-BN agrees well with 
the experimental value of 0.14±G.03 eV for the diffusion 
of Co adatoms on h-BN deposited on Ni(lll).^^ 

Figure [5l^a) shows the horizontal trend in the binding 
energies and the diffusion activation barriers across the 
third row of the periodic table. Both graphene and h- 
BN show double-peaked features with maximum binding 
energies at Ti and Ni which have (P and d^ electronic 
configurations, respectively, while the very weak binding 
of Cr and Mn corresponds to the situation of half-filled 
(i-shells. Strong binding energies with magnitudes up 
to 2.08 eV (1.14 eV) for graphene (/i-BN) indicate the 
contribution from covalent binding involving d electrons. 
However, the fact that the atomic magnetic moments due 
to the partially filled d-shells are largely preserved upon 
binding [Fig. ^c)] shows that graphene and /i-BN act 
as weak ligand fields with respect to the metal adatoms. 
Large magnetic moments associated to half-filled d-shell 
elements correspond to large values of spin splitting of 
d-electron states and, thus, to their reduced participa- 
tion to the covalent binding. The other elements in the 
first (second) half of the transition metal series tend to 
give enhanced (reduced) magnetic moments upon bind- 
ing. To elucidate the origin of this behavior, we per- 




FIG. 2: (Color online) (a) Horizontal and (b) vertical chemical 
trends in the binding energies Eh and the diffusion activation 
barriers _Ea of metal adatoms on graphene and ft-BN. Insets 
show the concerned parts of the periodic table: K-Ga in (a) 
and Co-Ir, Ni-Pt, Cu-Au in (b). (c) Horizontal and (d) ver- 
tical chemical trends in the magnetic moments M and the 
atomic charges Qm of the adsorbed metal adatoms in their 
lowest energy configurations. Crosses refer to the magnetic 
moments of isolated atoms. Electronic configurations of the 
isolated atoms are given at the bottom of the figure for refer- 
ence. 



formed the Lowdin population analysis^ and found that 
the adsorption of metal adatoms leads to the partial pro- 
motion of electrons from the s- to the d-shell thus leading 
to the observed changes of the magnetic moments. 

Such trends depending on the partial d-shell fill- 
ing appear to be generic since a very similar behav- 
ior was found for metal-benzene molecular complexesSS- 
and even for strongly bound substitutional impurities in 
graphenei^"— The binding energies on /i-BN as well as 
the activation barriers on both monolayers follow exactly 
the same qualitative trend within the transition metal se- 
ries (Sc-Zn). Outside of this region (K, Ca and Ga) only 
graphene is able to bind adatoms strongly since the inter- 
actions are then governed by the sole electrostatic con- 
tribution resulting from charge transfer. This contrast is 
well illustrated by comparing the Lowdin atomic charges 



Qm of metal adatoms on the two monolayers [Fig.[2jc)]. 

Vertical trends in the periodic table are illustrated 
in Fig. Hfb) for the late transition and coinage metal 
adatoms. In general, the binding energies tend to in- 
crease when moving down the columns of the periodic 
table. The reductions can be associated with the changes 
of the lowest energy adsorption sites discussed above. 
The largest binding energy of 2.02 eV is found for Rh on 
graphene while the highest activation barrier of 0.77 eV 
corresponds to Pt on /i-BN. At variance, the diffusion of 
Pd on graphene and of Rh and Ir on /i-BN are character- 
ized by very low diffusion barriers (<0.05 eV) combined 
with high binding energies (>1 eV). Inert coinage met- 
als have very low binding energies and diffusion barrier, 
with minimum values for Ag. Similar trends were also ob- 
served for carbon on coinage metal surfaces. '^^ Even on 
graphene, Ag shows practically no charge transfer as this 
element is located at the point of crossover between elec- 
tron donating (Cu) and electron accepting (Au) behavior 
[Fig-Hid)]. The hole-doping of graphene which is other- 
wise difficult to achieve by chemical means was recently 
demonstrated in Au deposition experiments.''"^ Interest- 
ingly, we find that also strongly binding Ir adatoms act 
as electron acceptors. 

It is worth stressing that the presently used theory 
does not provide a correct description of weak van der 
Waals interactions. This might quantitatively affect the 
results, especially when the considered adatom does not 
give rise to strong covalent bonding and/or electrostatic 
interactions. 



IV. EFFECT OF METALLIC SUBSTRATE 

On the surfaces of forth and fifth row transition met- 
als, epitaxial graphene and /i-BN produce Moire patterns 
due to the mismatch between the lattice constants of the 
substrate and of the monolayer [cf. Fig. [SJa)]. The lo- 
cal shifts of the monolayer lattice with respect to the 
substrate lattice across the Moire unit cell result in long- 
wavelength modulations of the PES which drive the self- 
assembly of periodic arrays of nanoparticles. These mod- 
ulations are due to variations in the monolayer-substrate 
distance and to changes in the local electronic structure.— 
Below, we focus on the role of local electronic structure 
variations in the diffusion of individual metal adatoms on 
metal-supported graphene and h-BN. Collective effects 
due to the local rehybridization induced by few-atom 
clusters have been considered elsewherci^^^ In addition, 
the distance modulations alone can lead to only minor 
curvature-induced effect on the adatom binding energies 
since the estimated height variations are actually rather 
small (A/i<1.5 A*!^) compared to the size of the Moire 
unit cell (a~3 nm)i^ 

In order to understand the effect of the local electronic 
structure variations, we study the binding of Co adatoms 
on epitaxial graphene and ft,-BN supported by lattice- 
matched Ni(lll) and Cu(lll) substrates. These two 
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FIG. 3: (Color online) (a) Schematic illustration of the Moire unit cell with three binding regions indicated, (b) Schematic 
representation of the potential energy surface (solid line) for the diffusion of the adatom moving along the horizontal direction 
of the unit cell shown in (a), based on calculations for Moire domains I-III. The actual local binding energies (dashed line) and 
diffusion barriers (oscillation amplitudes) correspond to the case of the Co adatom on graphene supported by Ni(lll). The 
dotted line refers to the case of the Co adatom on free-standing graphene. (c) Binding energy Eh of the Co adatom versus the 
charge transfer per carbon atom Qc from the metal to graphene for the three Moire domains. 



cases model the regimes of strong and weak monolayer- 
substrate interactions, respectively.-^'^ We stress that un- 
der realistic conditions these metals do not produce Moire 
pattern due to the small lattice mismatch<^i2i"— The 
three principal Moire domains [I-III in Fig. [Sja)] are 
modeled by introducing an artificial lateral shift between 
the monolayer and the substrate. Further details of this 
procedure are given in Sec. |TT1 

The calculated binding energies and activation barri- 
ers for Co adatoms on graphene|Ni(lll) show apprecia- 
ble variation across the three Moire domains [cf. Table H]. 
The resulting PES [schematically depicted in Fig. Elb)] 
can be viewed as a modulation of the PES of free- 
standing graphene [Fig.[lja)] by the long- wavelength per- 
turbations resulting from the periodicity of the Moire su- 
perlattice. The Moire domain I, which is characterized by 
the shortest graphene-substrate distance (2.23 A)^^— 
shows the lowest binding energy while the opposite is 
true for region II showing the largest graphene-substrate 
distance (4.11 A). However, Moire domain III presents a 
short graphene-substrate distance (2.49 A) in combina- 
tion with a high binding energy [Table |I] , indicating that 
there is no clear correlation between these quantities. At 
variance, a clear correlation can be established between 
the binding energies of the Co adatoms and the local 
charge transfer per carbon atom (Qc) from the metal to 
graphene [cf. Fig.lJfc)]. This result can be intuitively un- 
derstood: the charge transfer from the metal to graphene 
reduces the charge transfer from the Co adatom to the 
template surface and, thus, the electrostatic component 
of the binding energy. Very recently a similar electro- 
static effect has also been found for adatoms on ultra- 



thin oxide films. ■*° For Co adatoms on ft,-BN|Ni(lll), we 
found the same tendency, although the local variations of 
the binding energy are smaller in this case. On a weakly 
binding substrate such as the Cu(lll) surface, the PES 
remains practically unchanged across the Moire unit cell. 
The Moire domains associated with higher binding en- 
ergies are expected to show larger thermal populations 
of adatoms, and will thus act as nucleation centers for 
nanoparticles. Therefore, larger variations of the bind- 
ing energies across the Moire template will favor more 
ordered nanoparticle arrays with narrower size distribu- 
tions. We conclude that the templates made of graphene 
in combination with strongly binding metal surfaces (e.g. 
Ir, Rh, Ru) are more promising for the self-assembly of 
metal nanoparticle arrays. 



TABLE I: Comparison of binding energies Eh (in eV) and 
diffusion activation barriers _Ea (in eV) of the Co adatom on 
free-standing graphene (gr.) and /i-BN with corresponding 
epitaxial monolayers deposited on Ni(lll) and Cu(lll). The 
different positions of the monolayer atoms refer to three dis- 
tinct local regions of the Moire pattern shown in Fig. [S] 
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V. CONCLUSIONS 

We systematically studied the potential energy sur- 
faces of metal adatoms on graphene and /i-BN across 
the periodic table, and clarified the role of the metallic 
substrate in adatom diffusion. The present results are 
well understood within a simple picture involving cova- 
lent and electrostatic interactions. Our work formulates 
general principles required for the rational design of self- 
assembly templates based on epitaxial graphene and h- 



BN which can find applications in ultra-high density in- 
formation storage, catalysis and sensing, and so on. 
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